Some generalizations of the Raychaudhuri equation 
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The Raychaudhuri equation has seen extensive use in general relativity, most notably in the de- 
velopment of various singularity theorems. In this rather technical article we shall generalize the 
Raychaudhuri equation in several ways. First an improved version of the standard timelike Ray- 
chaudhuri equation is developed, where several key terms are lumped together as a divergence. This 
already has a number of interesting applications, both within the ADM formalism and elsewhere. 
Second, a spacelike version of the Raychaudhuri equation is briefly discussed. Third, a version of 
the Raychaudhuri equation is developed that does not depend on the use of normalized congruences. 
This leads to useful formulae for the "diagonal" part of the Ricci tensor. Fourth, a "two vector" 
version of the Raychaudhuri equation is developed that uses two congruences to effectively extract 
"off diagonal" information concerning the Ricci tensor. 

Keywords: Raychaudhuri equation, timelike congruence, spacelike congruence, non-normalized con- 
gruence. 
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The Raychaudhuri equation has become one of the 
standard workhorses of general relativity, particularly as 
applied to the singularity theorems. For textbook pre- 
sentations see for instance [11-13 ■ Nevertheless, we feel 
that there are still some interesting ways in which the 
general formalism can be extended. There are four ex- 
tended versions of the Raychaudhuri equation we wish to 
explore in this article: 

• Single timelike unit vector field. 

By collecting several terms in the usual formulation 
into a divergence, we obtain a particularly useful 
version that finds many applications in the ADM 
formalism and other situations. 

• Single .spacelike unit vector field. 

This situation is most typically ignored. We will 
make a few hopefully clarifying comments. 

• Single non-normalized vector field. 

This somewhat simplifies the Raychaudhuri equa- 
tion, at the cost of no longer having nice positivity 
properties. 

• Two non-normalized vector field.s. 

This allows us to probe off-diagonal components of 
the Ricci tensor. 

These four extensions of the Raychaudhuri equation will 
soon be seen to each be useful in their own way, and to 
yield diflFerent information. We shall provide numerous 
examples below. 
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This is the standard case. Let m° be a field of unit 
timelike vectors (a congruence). This does not have to 
be the 4- velocity of a physical fluid (though it might be), 
it applies just as well to the 4- velocities of an imaginary 
collection of "fiducial observers" [FIDOs]. Then it is a 
purelygeometrical result (see for example Hawking and 
Ellis "if, pp 82-84, or Wald 0, or Carroll Q, or Pois- 
son [4[, or even Wikipedia, (note that there are sometimes 
minor disagreements of notation — typically just a factor 
of 2 in odd places) that: 



as 



This is the standard form of the Raychaudhuri equation. 
The spatial projection tensor is 



hab = gab + UaUb. 



(2) 



This projection tensor has signature {0, +1, +1, +1}. 
Various shear and expansion related quantities are 



9 ^ab 



V„u": 



fafc — Qab — -^habO; 



^2 „ ab ^ r\ 

Vorticity related quantities are 



= WabW"'' > 



(3) 
(4) 

(5) 

(6) 

(7) 
(8) 
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With these definitions we have the usual decomposition 

1 dua 

Ua-b = i^ab + (yab + 7;^hab 7— Ub- (9) 

3 ds 

See (for example) pages 82-84 of Hawking and Ellis fH . 
Equation ^ is Wald's equation (9.2.11) [2], supple- 
mented with the Va (^n") term due to allowing a non- 
geodesic congruence; you can deduce the presence of this 
term from the second line in his (9.2.10) by not assuming 
geodesic motion. 

Now consider the identity 



u-V9 = V ■ [Ou) - 6'V • M = V • [Ou) 



(10) 



Using this identity we can also write the Raychaudhuri 
equation in the slightly unusual forms 



-Rabu'^u^' + uj^ -a^ + -e^, (11) 



RabU^u" = - ^2 + + (^-6 + . (12) 

This extension of the usual Raychaudhuri equation is 
"close" to, but significantly more general than, a key 
technical result used by Padmanabhan and Patel in ref- 
erences 0-01 ■ 



III. APPLICATIONS 

We now consider several applications of the above for- 
malism — these applications basically amount to strate- 
gically choosing an appropriate congruence. 



III.l. Vorticity-free congruence 

Let ^'(x) be an arbitrary scalar field and define a set 
of fiducial observers [FIDOs! by 



Ua OC Va^". 

Then normalizing we have 



|V*| 



and furthermore 



OJab = 0. 



(13) 



(14) 



(15) 



The minus sign here is purely conventional, it guaran- 
tees that the u° is "future-pointing" in the direction of 
increasing ^I'. Conversely 

UJab = =^ Ma OC V ■ (16) 

This is guaranteed by the Frobenius theorem. 



Then in this vorticity-free situation the extended Ray- 
chaudhuri equation reduces to 



as 3 \ ds 

or equivalently 



(17) 



Rabu'^u" = -a^ + -9^ + Va f -9u'' + ^ ) ' (18) 



or even 



du° 



RabU^u" = -eabO""' + e'+ Va \^-9u'' + ^ j . (19) 

But since u'^ is now hyper-surface orthogonal we can use 
the slices of constant ^E* to define a spacelike foliation — 
the scalar 4' serves (at least locally) as a "cosmic time" 
function. Then in terms of the extrinsic curvature Kab 
of the constant hyper-surfaces we have, (using Misner, 
Thorne, and Wheeler Q sign conventions for the extrinsic 
curvature), the results: 



9ab = -Kab-, 9 = -K; 



<Jab — - \ Kab — -^Khab 



" =2 



(20) 
(21) 

(22) 



But then 



Rabu^u" = -KabK-"' + K' + Va [K + —j . (23) 

This is effectively one of the key technical results used by 
Padmanabhan and Patel in [a-Q, but now we see that 
this result is actually a special case of a considerably 
more general result, and can be viewed as a relatively 
straightforward extension and then specialization of the 
Raychaudhuri equation. 



III.2. ADM formalism 

By definition, in any stably causal spacetime there is a 
globally defined "cosmic time" function t{x) such that dt 
is always timelike. Then on the one hand the constant-i 
slices are always spacelike and can be used to set up an 
ADM decomposition of the metric, while on the other 
hand u — —{dty/\\dt\\ is a vorticity-free unit timelike 
congruence, so that the results of the previous subsec- 
tion apply. (As usual, dt^ denotes the vector obtained 
form the one-form dt by "raising the index" , similarly 
will denote the one-form obtained from the vector u by 
"lowering the index.) 

Consequently the extended Raychaudhuri equation 
can now be cast in the form 



Rii = -KabK'^'' + 



^a[Ku'^+'^Y (24) 
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This result complements and reinforces the information 
one obtains from the Gauss equations — see for example 
Misner, Thorne, and Wheeler Q pp 505-520, or Ken- 
dall fo*] pp 23-24. The Gauss equations (for a spacelike 
hypersurface) are 

^^^Rabcd = '■''^Rabcd + KacKhd - KadKbc- (25) 

Contracting once 

= ^^'>Rab - ^^^RacMU^u" + tr(i^)i^„, - {K^)ab. 

(26) 

Contracting a second time 



= _ 2 '^^'^Rabu'^u' + K^- tr(X2). 



(27) 



But now, since ^^"^ Rabu"'u'^ has been given to us via the 
extended Raychaudhuri equation, we easily see that for 
a spacelike hypersurface 

Wi?= (3)i? + tr(x2)-/f2_2Va (^"" + ^) • (28) 

Traditional derivations of this result are sometimes some- 
what less transparent, and viewing it as an extension 
of the timelike Raychaudhuri equation is the cleanest 
derivation we have been able to develop. To see some 
of the deeper connections with the ADM formalism read 
(for example) §21.6 on pp 519-520 of Misner, Thorne, 
and Wheeler Q; note especially eq (21.88). See also ex- 
ercise (21.10) on p 519. Also note the discussion by by 
Padmanabhan and Patel in references [l-Q. Also, we 
should warn the reader that Wald uses an opposite sign 
convention for the extrinsic curvature. See specifically 
Wald equation (10.2.13) on page 256. 



III. 3. Static spacetimes 

Let us now take the discussion in a rather different di- 
rection, and assume that the spacetime is static. That 
is, there exists a hypersurface-orthogonal Killing vector 
/c" that is timelike at spatial infinity. Because it is hyper- 
surface orthogonal then ka oc Va^, and so — k°'/\\k\\ 
is a set of FIDOs of the type considered in the previous 
section. But since k"' is also a Killing vector we have 



k 



(a;b) 



and so obtain the quite standard result that 



^(a;6) 



V(.{A:/||fc||} 



k 



k{b^a) 



(a-b) 



(&Va)| 



(29) 



vorticity free) . The Raychaudhuri equation then special- 
izes to the particularly simple result 



R,, u-u' = V, ( ^ 
as 



(31) 



This is essentially the technical result we used in our 
derivation of an entropy bound for static spacetimes [lol . 
[Tl| . though in those articles we had derived it from an 
old result due to Landau and Lifshitz J,2] . (The original 
Landau-Lifshitz result is obtained via a straightforward 
but tedious series of index manipulations, with little ge- 
ometrical insight.) 



III. 4. Stationary spacetime — Killing congruence 

What can we now do for stationary, as opposed to 
static spacetimes? (This distinction is relevant to "ro- 
tating spacetimes" , for example Kerr spacetimes versus 
Schwarzschild spacetimes. See for instance [Tsl - flGj .) The 
(asymptotically) timelike Killing vector k = dt [that is, 
= (1; 0,0,0)"] is no longer hypersurface orthogonal. 
Nevertheless we can still define the timelike Killing con- 
gruence 



k" 

iM 



(32) 



This timelike congruence corresponds to a class of FIDOs 
[not ZAMOs, not zero angular momentum observers] that 
sit at fixed spatial coordinate position [l3, This 
timelike congruence, even though it is not hypersurface 
orthogonal, still satisfies equation (|29|) . So even though 
there is no longer any interpretation of the shear in terms 
of an extrinsic curvature, we still have 



whence both 



0, 



aab = 0; and 



Therefore 



RabU^'u'^u'+^a ( ^ 



(33) 



(34) 



(35) 



However, unless further assumptions are made, we can- 
not do much with the term. Generically we have 



U[a;b] = V[a{/c/||fc||}6] : 



k 



k\hS/„A\k\ 



Hence 



[fell ||fc||2 
'[a-b] W[fc||fc||,a] 



(36) 



7ab 



= 



Kab = 



K = 0. 



(30) This implies 



That is, in static spacetimes the extrinsic curvature of the 
time-slices is zero (in addition to the congruence being 



ab j^acj^bd 



(37) 
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whence 



RabU^U 



[a-b] K[c-d] 



\k\ 



(38) 



Unfortunately this does not simphfy any further, and 
without further assumptions for the timehke Kilhng con- 
gruence on a stationary spacetinie we should just be sat- 
isfied by the inequality: 



RabU^U^ > Va 



"d7 



(39) 



III. 5. Stationary axisymmetric spacetimes 

In a stationary axisymmetric spacetime consider the 
vorticity-free congruence of section IIII.ll (not the Killing 
congruence of section lTlI.4[) . Because of the axisymmetry 
the congruence u — —{dt)'^/\\dt\\ must then be a linear 
combination of the two Killing vectors, kt = dt and k^ = 
9^, in which case 9 — V ■ u = Q. In this case equation 
(ITSll reduces to 



Rabu'^u'' = -a^ + V J — 



du" 
d7 



(40) 



which implies, for the natural vorticity-free congruence 
on an stationary axisymmetric spacetime 



Rabu'^u' < Va 



dM° 
"d7 



(41) 



It is this particular inequality that we used in refer- 
ence [lH to place an entropy bound on rotating fluid 
blobs. (Note that the direction of the inequality has 
changed between equations (pQ)) and (|4T|) . but that is 
merely due to the fact that we are using different time- 
like congruences.) 



IV. SINGLE UNIT SPACELIKE VECTOR FIELD 

In counterpoint, we now let u° be a field of unit space- 
like vectors. The projection tensor becomes 



hab = gab - UaUb- 



(42) 



In contrast to the timelike situation the projection tensor 
is now of indefinite signature { — 1,-1-1,-1-1,0}. One can 
still formally define the quantities Oab, crab, and uJab, 
but they no longer have the same physical interpretation 
in terms of shear and vorticity. Furthermore since the 
projection tensor has indefinite signature we now can- 
not guarantee either cr^ > or cj^ > 0. On the other 
hand, the Raychaudhuri equation is formally unaffected. 
That is, the fundamental equations ([T]), ([TT|). and (HU, 
continue to hold as they stand. 

If we now consider a vorticity-free spacelike congru- 
ence, it will be hypersurface orthogonal to a timelike hy- 
persurface. (That is, the normal to the hypersurface is 



spacelike, while the tangent space to the hypersurface 
can be chosen to have a basis of one timelike and two 
spacelike tangent vectors.) 

In this situation we can without loss of generality set 
u = (d*)Vl|d*||. Then iOab 0, while in terms of 
the extrinsic curvature aab ^ ^Kab as for vorticity free 
timelike congruencies. Thus equation (|23p is formally 
unaffected and can now be cast in the form 



Rfi 



-KabK-" 



du" 

"d7 



(43) 



On the other hand, because u is now a spacelike normal 
to a timelike hypersurface there is a key sign flip in the 
Gauss equations, which now read 



^^^Rabcd — ^^^Rabcd — KacKbd + KadKbc- 

Contracting twice 

Therefore for a timelike hypersurface we have 



du'' 
"d7 



(44) 



(45) 



(46) 



In summary, for spacelike congruences the Raychaud- 
huri equation itself is formally unaffected (though the 
projection tensor is slightly different and we can no longer 
rely on the non- negativity of and o;^). However appli- 
cations of the Raychaudhuri equation, specifically any- 
thing involving the Gauss equations for embedded hy- 
persurfaces, typically exhibit a limited number of sign 
flips. 



V. SINGLE NON-NORMALIZED VECTOR 
FIELD 

Now consider an non-normalized vector field u° , either 
spacelike, timelike, or null. What if anything can we say 
about the quantity 

??? (47) 

Following and modifying the discussion of Wald , see 
(E.2.28) on page 464: 

RabU'^u'' = R^acbu'^u'' 

= -U" [VaVfc - VfcVa] 

= ^Vaiu'^Vbu'') + (VaU°)(Vfcu'') 

+ Vfc(u"Vau')- (Vfcu")(Vau'') 
= Va(-U°V6U^ + u''V6U°) 

+ (V-u)2-(Vfc«a)(W) 

= V • {(u- V)u- (V • u)u} 

+ (V-li)2-(V6Ua)(W). (48) 

In obvious notation, using 6 — \7 ■ u, this can be cast as 



RabU^u" = V • {V„w - 0u} 



(49) 
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This result can be viewed as another generalization of the 
Raychaudhuri equation. The advantage of this particular 
formula is that we have not carried out any projections, 
and have not even committed ourselves to the nature of 
the congruence, be it spacelike, timelike, or null. One 
disadvantage is that because of the Lorentzian signature 
of spacetime we cannot (at least not without further as- 
sumption) guarantee 



??? 



\/^aUb] vI'^m''! > 111 
Two specific applications come readily to mind; 



(50) 



(51) 



• For any Killing vector u°- — fc" we have V(QMf,) = 0, 
and consequently 9 = Q. Therefore for any Killing 
vector whatsoever 



Rabk^'k'' = V • {Vfefc} + V[,A:fc]Vl"fc''l. 



(52) 



• For any one arbitrary exact one-form u = c?^', even 
a locally exact one-form, we have V^aUftj — 0, while 

e = V^^- and V(aU6)V("M^) = ^';a;&^''°'^ There- 
fore for any locally exact one-form whatsoever 



i?°''(d*)a(d^')h = V • {Vdfd* - (V^ 4')d«'} 



(53) 



In fact, ^ could simply be one of the spacetime 
coordinates (defined on some suitable local coordi- 
nate patch) in which case this version of the Ray- 
chaudhuri equation turns into a statement about 
the diagonal components of the Ricci tensor in a 
coordinate basis 

i?** = V • {Vd^d^- - (V^ «')d*} 

+ (V2il»)2 _ (54) 

More boldly, if one chooses ^' to be a harmonic 
coordinate, (V^^ — 0), and this can always be done 
locally, then 



(55) 



In summary, this extension of the Raychaudhuri equation 
has given us some useful computational formulae. 



VI. TWO NON-NORMALIZED VECTOR 
FIELDS 

We shall now ask if it is possible to extract any use- 
ful information by considering two different congruences 
simultaneously. 



VI. 1. Motivation 



To motivate this particular extension of the Raychaud- 
huri equation, recall that many decades ago Landau and 
Lifshitz had shown that in any stationary spacetime [l3| 
(§105 equation (105.22), for a recent application of this 
result see [13, Ell): 



-54 



Oa T^-i \ 

-gig r ao) 



(56) 



(Here a € {0,1,2,3}; i G {1,2,3}.) But because the 
metric is stationary {t independent) we can also write 
this as 



1 



-54 



(57) 



To begin converting this into a coordinate-free statement, 
note that 



i?o" = R\{dt)a{dtf = R\ {dt)a k" 



(58) 



Here we have had to use both the timelike Killing vector 
k, for which fc"^ = (9t)" = (1,0,0,0)", and the one-form 
dt, for which {dt)a — (1,0,0, 0)a. By direct computation 



5°° r^o = g'" r\d {dt),k'' = r^^ {dtyk'' 

= r^rf k^dtf = {ddk" + T\dk'} [dty 



(59) 



But then 



Ro 



-.94 



-545 r aO) 



1 



db (y^(d<)'*(Vdfc'')) 



V-.94 
VbiidtfiVdk')}. 



(60) 



So the Landau-Lifshitz result is equivalent to the state- 
ment that in any stationary spacetime 

R\ (dt)afc^ = WbiidtyiWdk')} = V • (Vdt.fc). (61) 

So some linear combination of Ricci tensor components is 
given by a pure divergence. Note that two different vec- 
tor fields are involved. This observation naturally leads 
to the question: Is it possible to come up with a vari- 
ant of the Raychaudhuri equation that depends on two 
congruences and u"? Something of the form 



Rab m'w' 



??? 



(62) 



We shall see how this is done below. 
For now, let us mention that 

{Vdk'Kdtf = (v^fc'')(dt)d = -{v'k''){dt)d 

= -V''{k\dt)d} + k''v\dt)d 
= -V''{1} + fc'^V^Vdi 
= k'^'^^Vdt = k-^VdVh. 



(63) 
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So the Landau-Lifshitz result can also be written in the 
alternative form 

R\ (di)afc' - Vbik'^VdVh} - V • (Vfedt«). (64) 

Finally note that 

(ydk''){dt)\dt)t, = {^dkb){dtY{dtf = 0, (65) 

so the vector V^tjA: ~ Vfcdt' is perpendicular to dtK 

VI. 2. Construction 

Following and modifying the discussion of Wald , see 
equation (E.2.28) on page 464: 

= -U" [VaVh - VhVa] 

= -Vaiu'^^bv'') + (V,M'^)(VfcW^) 

+Vb(u"Va«'')-(Vbii")(V„«^) 
= Va(-?i°V6w'' + w^Vbi;") 

+(V-7i)(V-i;)-(Vbii,)(VV). (66) 

With minor notational changes and given the symmetry 
of the Ricci tensor this can also be written as 

RabU^v^ = V • {(m • V)v - (V • v)u} 

+ (V-u)(V-«)-(V6u,)(V°f''), (67) 

and 

R^bU^v^ = V • {(w • V)u - (V • u)v} 

+ {V ■u){V ■v)-{ybUa){V''v''). (68) 

Furthermore (in obvious notation) this can again be 
rewritten as 

Rabu^v^ ^ V • {V„w - e^u} + OuOv 

-V^aUb)V^^v''^ + V[aUb]V^^v''l (69) 

and 

Rabu'^v'' = V • {V„M - OuV} + OuOv 

-V(ali6)V(^7'')+V[„U6]V[^7''l. (70) 

Note the similarities to the single-congruence case, and 
note particularly the presence of a divergence term gen- 
eralizing the standard Raychaudhuri equation. To check 
the equivalence of these two formulae note 

[VuV - Oyu) - (Vi,u - 9uv) 

= i^uV + Ouv] — [VvU + 9yu\ 

= V- [u(^v — v(^u] 

= \/-[uAv]. (71) 

That is, the difference of these two currents is the diver- 
gence of a 2-form, which makes it automatically closed. 



VI. 3. Generalizing the Landau-Lifshitz result 

Let u = fc be any Killing vector, and let be an 
arbitrary (locally) exact one-form, so w = (dvP)" where 
\l/(x) is an arbitrary scalar. Then V(oWf,) = '^(akt) = 0, 
and so f?„ = 0. Furthermore ViqU;,] = V[aVh]4' = 0, so 
from equation (j69p we have 

fc^'V^* = V • {Vfcd* - (V2*)fc}, (72) 

while from equation ()70p we have 

i?a6fc"V^*- V-{Vd*«A:}. (73) 

This nicely generalizes the Landau-Lifshitz result to any 
arbitrary Killing vector and any arbitrary (locally) exact 
one form d^, not just dt. (That these two formulae are 
equivalent follows from the discussion in the previous sec- 
tion above.) Note the (standard) Landau-Lifshitz result 
corresponds to — (9*)" and — 

Now choose a coordinate system adapted to the Killing 
vector k. Let k = 9k define a Killing coordinate K, so 
that all geometrical objects are independent of the coor- 
dinate K. Let also be viewed as a coordinate, relabel 
it as a;", possibly distinct from K, and with no claim that 
necessarily corresponds to a Killing vector. Then 

i?K" = V-{V(d,a).9K}. (74) 
Unwrapping the covariant derivatives we see 

Rk"" = -iL= db (V^s'T^k) . (75) 

If we now let the index i range over every coordinate 
except the Killing coordinate K then, because all geo- 
metrical objects are independent of the coordinate K, we 
have 

i?K" = a,; (\/^5°' T^k) . (76) 

This equation, ultimately based on our two-congurence 
extension of the Raychaudhuri equation (1691) , is now very 
much in Landau-Lifshitz form, but us definitely consid- 
erably more powerful than the original Landau-Lifshitz 
result. 



VI. 4. Landau-Lifshitz in cixial symmetry 

Since in a stationary spacetime with axial symmetry 
we have a second azimuthal Killing vector — >■ (d^)"", 
and could also consider — ^ 0, then there are three 
additional Landau-Lifshitz like results: 

i?/ = i?afc(a0)°V^i = V-{Vdt.94; (77) 
i?/ = i?,b(9t)"VV = V-{Vd^.aj; (78) 
i?/ = i?ab(9^)"V''(/) = V-{Vd0«90}. (79) 
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Let the indices A, B G {t, (/>} then we can cohect these 
results (four of them altogether) as 



Unwrapping the covariant derivatives 
1 



-54 



db{V^4g'"'r''aA) 



(80) 



(81) 



If we now let the index i range over every coordinate 
except the two Killing coordinates t and (f>, then 



Ra^^ 



1 



'-94 



(82) 



Making this all very explicit, there are now four Landau- 
Lifshitz like results in total. They are: 



Rt'^ 



1 



-54 



-54 5*"r\0 



(83) 



diagonalize the metric: gab = 9ab ffi 9ij- But then 
g""" T\a = g"""" T'cA = g''''g''V,cA 
= -Ig"""" 9'' d.gcA- 
So finally we have the relatively compact result 
1 1 



(87) 



Ra" = - 



2 



^^{^4g'''' g'' d.gcA) 



This can be rearranged in a number of different ways. As 
an illustration we point out 



Rab 



1 1 



2 



'94 g'^ djgAB) 



+ 2 5'^ ^^9AC 9 djgoB- 



(89) 



We again see that our two-congruence extension of the 
Raychaudhuri equation has given us additional useful in- 
formation regarding the Ricci tensor which might be dif- 
ficult to extract by other means. 



Rt'^ = -1= ^^ (V^5^" ^\t) ; (84) 



VII. DISCUSSION AND CONCLUSIONS 



In this somewhat technical article we have developed 
several useful extensions of the usual Raychaudhuri equa- 
^ tion. The main theme has been to relate various linear 

— r)i 54 5*° r*o0) ; (85) combinations of components of the Ricci tensor to di- 

V 94 vergences of suitably defined fluxes. We have worked 

with timelike congruences, spacelike congruences, and 
non-normalized congruences, in all cases being able to 
(p _ ^ g f ■\/ ~54 5*^° ) (86) j^^^ ^ little bit more (and sometimes much more) 

\J—g4 than standard the Raychaudhuri equation would imply. 

One potentially far-reaching result is the "two congru- 
Furthermore, recall that in stationary axisymmetric ence" extension of the Raychauduri equation presented 
spacetimes we can always choose coordinates to block in equations and (fTO]) . 
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